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""Ji; ABSTRACT 

, The star-to-star differences in the abundance of light elements observed in the globular 

■ clusters (GCs) can be explained assuming that a second generation (SG) of stars form 

in the gas ejected by the asymptotic giant branch (AGB) stars belonging to a first 
stellar generation. However, while Na and O appear to be anticorrelated in the cluster 
stars, from the stellar models they turn out to be correlated into the AGB ejecta. 
Q I In order to reconcile the stellar theory with the observational findings, all the GC 

models invoke an early dilution of AGB ejecta with pristine gas occurring during the 
C/3 ' SG formation. Despite a vast consensus about the occurrence of such dilution, the 

^ , , physical process behind it is still unknown. In the present paper we set some general 

constraints on the pristine gas dynamics and on the possible amount of pristine gas 
' involved in the SG formation, making use of a one zone chemical model. We find 

^ , that such a dilution is a necessary ingredient in the SG star formation to explain 

lO ' the observed abundance patterns. We confirm the conclusion of our previous works 

showing that clusters must have been initially much more massive. We also show 
. that models assuming that clusters had an initial mass similar to their current one, 

' and adopting a large fraction of pristine gas to form SG stars, fail to reproduce the 

ff^ I observed Na-0 anticorrelation and are not viable. We finally show that the dilution 

i event should be restricted in time, rather than extended for the all duration of the SG 

formation. 
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1 INTRODUCTION currently on the main sequence (MS). This implies that the 

above chemical characteristics have been imprinted in the 
In the last decade a number 01 spectroscopic data provided , . r ^ i i mo 

gas by a previous generation oi stars, because low-mass Mb 
the evidence of spreads of hght elements and anticorrela- , , . ^ iii 

° r , n stars do not reach the high temperatures needed to sustain 

tions between Na and O and Mg and Al m the stars or glob- ^, , ^- t t ^ ^^ , i i i 

r n the nuclear reactions leading to the observed relations be- 

ular clusters (GCs), indicating the presence oi at least two , i ^ 

^ ' . tween the elements, 

stellar populations in these objects. As these properties are 

not shared by the halo stars with the same metallicity, nor 

by open clusters in our or other g alaxies, they represen t a it is still debated which kind of stars act as pol- 
trait of the GCs so peculiar to leadlCarretta et all l|2010l) to Inters, whether asymptoti c giant branch (AGB) stars 
identify GCs as those clusters where there is a Na-O anticor- (|DAntona fc Ventura"2007') or fast rotating massive stars 
relatio n. On the basis of this anticorrelation, [CajrettaelalJ (FRMsTt'ecressin ct al. 200lB). In|p'Ercolc ct al. (200^), 
1)2010') separate the GC stars in three families: a Primordial we presented a model for the formation and dynamical evo- 
population of stars sharing the same abundances of the halo lution of multiple populations in GCs in which second gen- 
stars; an Intermediate population of stars with a lower O eration (SG) stars form out of the AGB ejecta of the first 
abundance and higher Na; and an Extreme population of generation (FG) stars. We first focussed on the hydrody- 
very O-poor stars ([0/Fe]< -0.4). namic of the ejecta that collect into the cluster core through 
The variations are also found among unevolved stars a cooling flow; it turned out that, in order to shed enough 

gas to form the observed SG, the FG should have been orig- 
inally about ten times larger than that observed today. By 

* E-mail: annibale.dercole@oabo.inaf.it means of N-body simulations we then showed how the clus- 
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ter can lose FG stars reaching the nowadays observed ratio 

of the number of FG and S G stars. 

In a successive paper (jD'Ercole et all l2010l . hereafter 
Paper I), we focused on the chemical evolution of the GCs 
through simple one-zone models, and we were able to repro- 
duce the main characteristics of clusters such as NGC 2808 
and M4. The main problem encountered in this kind of mod- 
els is given by the fact that the current stellar models indi- 
cate that during the AGB phas e [0/Fe] and [Na/Fe] in the 
e.iecta are always correlated (e.g. Karakas fc Lattanzidl2007l: 



Ventura fc D'Antonal2009l : ISiessll2010l : IVentura fc D'Antonal 



2OIID . Thus, the O and Na abundances in the cluster inter- 
stellar medium (ISM) are correlated, and the SG stars form- 
ing in this gas would exhibit the same correlation, in glaring 
conflict with the observations. Also in the FRMS scenario 
some degree of dilution with pristine gas is needed to account 
for the presence of lithium in SG stars bec ause this fragile 
element is absent in the massive star eject a (|Decressin et al.l 
l2007al : iLiiid et al.llioTll '). To date, dilution seems to be re- 
quired by all the proposed models. 

Dilution was actually implemented in the models pre- 
sented in Paper I. This implementation, however, was real- 
ized via a reasonable but arbitrary modelling. In fact, de- 
spite the general consensus about the need of dilution, the 
source and dynamics of the diluting pristine gas is still not 
clear, and the mechanisms proposed t ill now ( D'Ercole et al. 
20081 : iBekki fc Mack cv 2009; Pflamm-Altenburg fc Kro upa 



20091 : IConroY fc Spergellboiol : iGratton fc Carrettall2010l) aU 



suffer from some drawbacks. 

The aim of the present paper is to shed some light on 
the dilution process through possible constraints that can 
be inferred by the chemical properties of the GCs. 



2 ON THE NEED OF DILUTION 

Given the vagueness of the dilution process and the great 
uncertainties plaguing some relevan t aspects of the nu- 
clear p hysics during the AGB stage jVentura fc D'Antonal 
l2005al lbl. Paper I) ,one could be tempted to make the simple 
(but audacious) statement that the actual stellar models are 
incorrect, and that the Na abundance in the AGB atmo- 
sphere decreases in pace with the decrease of stellar mass. 
In the end, the O-Na anticorrelation shown by the SG stars 
would be the direct consequence of the O-Na anticorrelation 
already present in the gas in which these stars form. 

Unfortunately, although the above picture eliminates 
some difflculties connected with the dilution, other problems 
remain unsolved. While the O-Na anticorrelation seems to 
be a distinctive characteristic of GCs, the e xtension of this 
antic orrelation varies from cluster to cluster l|Carretta et al.l 
I2OO9I ). In the previous scenario the variation of such exten- 
sion - that is, of the minimum value [0/Fe]inin El - depends 
on the time at which the star formation of the SG starts. 
For the sake of example, we can examine a "fake" evo- 
lution of [Na/Fe] within the gas shed by the AGB stars, 
with the Na abundance decreasing in time and the oxy- 
gen which instead grows. Let us now consider two clusters. 



GC#1 and GC#2, starting to form SG stars at the times 
t\ and t2 > ti, respectively. Clearly, GC#1 will host some 
very O-poor stars which instead will be absent in GC#2, 
and therefore the anticorrelation in this latter cluster will 
be shorter. On the other hand, given our assumption about 
the Na evolution, the maximum value [Na/Fe]max will be 
larger in GC#1 than in GC#2. Thus, an anticorrelation 
is expected between [0/Fe]min and[Na/Fe]niax, at odd with 
th e observations showi ng instead a correlation (see Fig. 19 
in lCarretta et al.llioogi '). 

While this correlation is not quantitatively accounted 
for by any model yet, we think that in any case it falsi- 
fies the starting hypothesis that the anticorrelation O-Na is 
imprinted in the stellar models. Therefore, despite the prob- 
lems outlined above, we more safely acknowledge the current 
models of stellar evolution, and believe that the dilution is 
an inescapable ingredient of the GC evolution. 



3 CONSTRAINTS ON THE FG MASS 



As discussed in iD'Ercole et al.l l|2008l 'l. in order to ex- 
plain the current mass of SG stars, globular clusters must 
have been initi a lly ab out ten times larger than toda}0. In 
IVesperini et~all (|2010l ). we have further explored the actual 
amount of SG stars that a cluster could have formed depend- 
ing on its initial mass and structural parameters (along with 
the implications for the contribution of globular cluster FG 
and SG stars to the assembly o f the G alactic halo). 

However, iDecressin et al.l (|2007al ) leave open the pos- 
sibility, in the FRMS scenario, of an initial FG mass 
close to the present one, provided a ra ther flat IMF. The 
same hypothesis was put forward by iD'Antona fc Caloil 
assuming the AGB stars as polluters. As shown by 
ID'Ercole et"al] ((20081 '). a flatter IMF increases the amount 
of available ejecta to form the SG population at the cost of 
reducing the long-lived FG stars, and a large Mfg is still 
needed to raise their number; as a consequence, a flat IMF 
is not a viable hypothes is to avoid larger clusters in the past. 

More recently also IConrov fc Spergel (|2010l 'l have pro- 
posed a model, still based on the AGB scenario, in which the 
cluster has an initial mass similar to the current one; the SG 
stars form from the small amount of AGB ejecta available 
and a very large amount of pristine gas accreted by the GC 
while moving through the ambient ISM . This last scenario 
can be easily ruled out for those few clusters in which a very 
high helium enrichment {Y ~ 0.38) has been safely derived 
from their CMDs because such a large amount of helium 
can not be released by a number of FG stars as small as the 
present one ([Renzini 200S|). 

Here we show that, apart from the helium abundance, 
the assumption of a FG comparable in mass to the present 
one looks impractical on the basis of the observed anticor- 
relation between O and Na. Let us assume, as in Paper I, 
a SG star formation starting at t = 32 Myr and lasting up 
to a fiducial time tend ~ 100 Myr. During thi s lapse of time 
the F G of mass AIpc with a Kroupa IMF (jKroupa et al.l 
Il993l ) sheds an amount of ejecta A/cj = 0.05Mfg (see 



^ The extension of the anticorrelation is regulated by the value 
of [0/Fe]min as the maximum value is essentially constant in all 
the clusters and is the same found in the halo stars. 



The actual factor depends on the stellar IMF, the duration of 
the SG for mation event, the am ount of ejecta a cluster is able to 
retain 



(see IVesperini et al.ll2010l) 
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iD'Ercole et al.|[2008l ). If we now further assume to form an 
amount of SG stars corn parable to the FG, as suggested by 
IConrov fc Soergell the quantity of the accreted pris- 

tine gas must be of the order of Mpr ~ Mpg- In conclusion, 
the degree of pollution of the SG is about Mej/Mpr = 0.05, 
which is definitively too low to explain the observed exten- 
sion of the O-Na anticorrelation 0. 

The above conclusion is highlighted by our chemical 
one-zone model shown in in Fig. [1] . Adopting the formalism 
of Paper I, the model is characterized by the following set of 
parameters:(tond,7, P*,fg, /pr, x) = (10, 700, 1, 1, 0.5). Here 
tend, 7 represents the time at which the evolution stops in 
units of 10^ yr, p*,fg is the FG density in Mq pc~'^, /pr is the 
amount of accreted gas in unit of /5*,pg, i' is the star forma- 
tion efficiency, and x is the ratio between SG and FG stars. 
Contrary to the models shown in Paper I, where the accre- 
tion rate of the pristine gas is assumed to have a gaussian 
temporal profile, here we set a c onstant accretion rate, as 
expected in the case envisag ed bv lConrov fc Soergell ()201GI 1 
of a GC moving with constant velocity through an uniform 
medium. 

The top panels of figure [T] show that the FG and SG 
stars occupy the same region in the [Na/Fe]-[0/Fe] diagram 
and share the same helium distribution extremely peaked 
around Y = 0.246, the value assumed for the pristine gas 
and the FG stars. The two populations are thus undistin- 
guishable. 

Inhomogeneous mixing has been suggested to overcome 
the above shortcoming and to ob tain the observed anticor- 
relations ijConrov fc Spergelllioiol ). However, a simple argu- 
ment shows that resorting to inhomogeneous mixing does 
not solve this problem. Let us consider an amount of AGB 
ejecta Mcj with a typical sodium abundance [Na/Fe]=0.8 
merging with a mass Mpr of pristine gas with, for exam- 
ple, [Na/Fe]=0. In case of inhomogeneous mixing, differ- 
ent regions of the mixture would have different values of 
[Na/Fe], and the stars forming in this gas would have a dis- 
tribution function A'^([Na/Fe]) extending across the range 
< [Na/Fe] < 0.8. Although the shape of this function de- 
pends on a number of assumptions (see the Appendix), some 
concl usions can be drawn i ndependently of its exact knowl- 
edge. ICarretta et al.l (|201G| ) separate the FG from the SG on 
the basis of their sodium abundance, with the SG stars hav- 
ing [Na/Fe] > 0.43. Thus, assuming that all the mass Afcj of 
the AGB ejecta forms a number A^'tot of stars, only a fraction 
^ = ([Na/Fe] > 0.4)/Ntot of them wiU appear as SG stars, 
the rest remaining confused with the genuine FG stars. The 
apparent ratio between first and second generation therefore 
will be: 

Msg _ AfsG,[Na/Fc]>0.4 _ CC /-.n 

Mfg ^ MpG + AfsG,[Na/Fc]<0.4 " 1 + (1 " C)C' 

In the above equation we have assumed that the mass lost 
by the FG stars (before t = 100 Myr) isAfoj = .^Mfg. From 



the observations we know that Msg /Mfg ^ 1, and therefore 
equation [T] gives rise to the condition 

0-5-!-y^ < C < 1- (2) 

This inequality is fulfilled only for ^ > 1, which is clearly 
meaningless q 

In summary, if the initial FG mass was similar to the 
present one, the observed O-Na anticorrelation cannot be 
reproduced in any way by inhomogeneous mixing. 

To give a "visual" support to the above conclusion, we 
have implemented an algorithm in our one- zone code in or- 
der to take into account inhomogeneous mixing (see the Ap- 
pendix). The results are shown in the lower panels of Fig. 
[1] Non homogeneous mixing clearly acts in the right direc- 
tion, but the number of Intermediate stars turns out to be 
rather small because of the large overabundance of pristine 
gas compared to the AGB eject4f|. This is particularly ap- 
preciable in the helium distribution, which acquire only a 
small tail toward higher values of Y , although the ejecta of 
massive AGB stars have values as high as V = 0.38. 

We thus conclude that a large FG mass must be a 
prerequisite of every model of GC evolution, which then 
must account for ho w to get rid of the stellar surplus (see 
iD'Ercole et"al]|2008f ). 

At the end of this section, we mention the effect of a 
possible inhomogeneous mixing in models similar to those 
of Paper I, with an initial Mfg much larger than today, 
and a; ~ 1. In this case it turns out that, in order to create 
the same amount of recognizable SG stars, inhomogeneous 
mixing models require a FG mass 10-70% larger than in the 
homogeneous case (see the Appendix). 



4 CONSTRAINTS ON THE TEMPORAL 
EVOLUTION OF DILUTION 

In the previous section we considered the case of a globular 
cluster moving in the halo through the ambie nt medium that 
can b e captured via the Bondi accretion (e.g. iLin fc Murravl 
'2007*). In addition to this process, the moving cluster may 
accrete further material sweeping up the surrounding gas 
which impinges on a seed of AGB ejecta present in the GC 
centre. Both the above accretion mechanisms may work only 
for a limited range of values of the parameters of interest. In 
fact, if the energy of the incoming gas is much larger than 
the GC potential well, a collective accretion is prevented. 
Moreover, if the cluster velocity is too large, the ram pres- 
sure exerted on the gas at the bottom of the potential well 
becomes stronger than the gravitational restoring force, and 
the sweeping process clears out the cluster of its gas instead 
of promoting accretion (see lLin fc Murravll2007^ . 

A full hydrodynamic study of the processes depicted 
above will be presented in a separate paper (D'Ercole et 
al. in preparation). Here we want to constrain the accretion 



^ Note that we have minimized the dilution posing Msg /A^fg = 
1; assuming a larger ratio MgQ / M-pG (^^ usually observed) would 
lead to an even more irrelevant pollution. 

^ Actually, the separating value is not the same for all the clus- 
ters, but is 0.4 dex larger than the minimum value of [Na/Fe] 
observed in the cluster. 



^ In models in which the FG is at least ten times more massive 
than the SG, the above condition becomes C, < 0.1(1 + £.)/£., and 
is easily met for any value of 

^ in order to maximize the effect of the inhomogeneity, we as- 
sumed a = 1.5, a being a parameter regulating the shape of 
Af([Na/Fe]) (see the Appendix). 
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Figure 1. Model with uniform accretion with the following 
parameters:(tcnd,7, P*,FG> /pr, I', a;) = (10,700,1,1,0.5) Top-left 
panel: chemical evolution of the homogeneous model in the plane 
[0/Fe]-[Na/Fe]. In this panel the dot-dashed line represents the 
ISM chemical pattern in absence of any accretion of pristine gas, 
and the dashed line (barely visible) the actual GC evolution for 
the assumed dilution. The FG stars are represented by the blue 
circles, and the SG stars by the cyan circles. Top-right panel: 
helium distribution of the homogeneous model. Bottom-left and 
bottom-right panels are similar to the upper ones, but in case of 
inhomogeneous mixing. 

process on the basis of chemical arguments. If GCs actually 
accrete gas from the ambien t medium continuo u sly sin ce the 
beginning, as suggested by IConrov fc Spergell \2Ql(h . dilu- 
tion starts quite soon and the formation of the Extreme pop- 
ulation is inhibited because the stars of this population are 
extremely 0-poor, and therefore are expected to be formed 
in the pure ejecta of the massive AGB stars. 

In order to explicate this point, we run a model simi- 
lar to that showrn in the previous section, but assuming an 
amount of accreted matter which is only 5% of the FG mass, 
a characteristic value in the models shown in Paper I: we 
thus assume (fend,?, P*,fg, /pr, J^, a;) ~ (10,700,0.05,1,0.5). 
As expected, in the homogeneous mixing case Fig. [2] (up- 
per panel) shows that the stars have all [O/Fe]>0. Once the 
inhomogeneous mixing is allowed (see Fig. [2] lower panel), 
some 0-poor stars actually form, but they remain too few 
to fit the amount of the Extreme (and Intermediate) stars 
in the G C where they are observ ed, reaching 15-25% of 
the total (jP'Antona fc Caloilliool ). This is well illustrated 
by the helium distribution function which barely reaches 
Y ~ 0.33 in the mo del, while nearly 20% of the stars in NGC 
2808 has Y > 0.35 (jCarretta et alJl2009l : lD'Antona fc Caloil 
|2008| ). Moreover, comparing the two distribution functions 
in Fig. [21 we note that in the inhomogeneous case the peak 
at y ~ 0.25 given by the FG stars is higher, reducing the 
ratio Msg/Mfg with respect to the homogeneous case. This 
is an effect due to the spread of the SG population which 
partially overlaps the Primordial stars (see the discussion in 
the Appendix). 



Figure 2. Model similar to that illustrated in figure [T] but with 
/pr = 0.05. Top panels: homogeneous case. Bottom panels: inho- 
mogeneous case. 

The above argument rules out not only continuous ac- 
cretion of external pristine gas, but also the internal c ontin- 
uous dilution suggested bv lGratton fc Carrettal (|2010l ) with 
the pristine gas produced by mass lost from less evolved stars 
(typically - but not only - MS stars) than those producing 
the polluting material. 

We thus conclude that the dilution can not be achieved 
by a continuous mixing of pristine gas, at least for those 
clusters showing the presence of an Extreme population. We 
remind here that the models shown in Paper I may form O- 
poor stars because the bulk of the accretion occurs at a 
time tac larger than the time at which the massive AGB 
stars start to pollute the cluster. Therefore, accretion seems 
to be a quite episodic event, rather than a regular one. 



5 CONCLUSIONS 

In this paper we highlighted some characteristics of the di- 
lution process occurring in the CGs on the basis of their 
chemical properties. Our findings can be summarized as fol- 
lows: 

(i) In principle, one could reproduce the ubiquitous O- 
Na anticorrelation present in the GCs without any dilution, 
simply assuming that all the available AGB models are in- 
correct and that actually the ejecta of the smaller AGB stars 
are less abundant in sodium and more abundant in oxygen. 
However, we have shown that this option appears to be ruled 
out by th e direct correla, t ion be tween [0/Fe]min-[Na/Fe]max 
found by ICarretta et al.l (|2009l ) . We thus conclude that di- 
lution must have actually occurred in GCs. 

(ii) We have studied the dependence of the extent of the 
Na-O anticorrelation on the amount of pristine gas involved 
in the SG formation. Our analysis shows that the contribu- 
tion of pristine gas to the matter from which SG stars form 
can not be larger than ~ O.IMfg. This confirms one of the 
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conclusions of our previous papers (Paper I; iD'Ercole et alj 
I2OO8I : FVesperini et al.ll2oTol ): in order to form the amount of 
SG stars observed today, clusters must have been initially 
about ten times more massive. Models suggesting an ini- 
tial mass similar to the present one, and based on a flat 
ISM (jP'Antona fc Caloill2004l: IPecressin et al.ll2007"al ). have 



been ruled out by D'Ercole et al.l (I2OO8I '). A model recently 
proposed bv lConrov fc Spergell (|2o"l?)" suggests that clusters 
would not need to be initially much more massive than to- 
day since SG stars would form almost entir ely from a larg e 
amount of accreted pristine gas. Following iRenzinil (|2008l) , 
such a model can be easily ruled out for those massive clus- 
ters containing a very helium rich population, because a clus- 
ter with a FG mass initially similar to the present one can 
not release the needed helium mas s. The results presented in 
this paper imply that the model of lConrov fc Spergell (|201Cll ) 
is ruled out also for all the other clusters without any Ex- 
treme Helium-rich population. In fact, if the accreted mass is 
of the same order of the FG, the abundances of the SG stars 
are almost identical to those of the FG stars, thus preventing 
the formation of the O-Na anticorrelation which instead is 
present in all GGs. Even a possible inhomogeneous mixing 
between the AGB eject a and the ac creted pristine gas, as in- 
voked bv lConrov fc Spergell (|2010l ). can not give rise to the 
Extreme SG population and to a substantial Intermediate 
SG population. 

(iii) We demonstrated that if the dilution starts at the 
same time of the pollution, the Extreme stars can not form, 
because such O-poor stars are presumably built in the pure 
ejecta of massive AGB stars. Thus, some proposed accre- 
tion processes such as th e Bondi accretion and streaming 
ijConrov fc Spergell |2010|'). or the self-dilution due to un- 
evolved stars ( Gratton fc Carrettall20l"ol '). can not work, at 
least for clusters hosting an Extreme population. Instead, di- 
lution should be considered as a delayed episode restricted 
in time, as described in Paper I. 
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APPENDIX A: INHOMOGENEOUS MIXING 

In order to account for a possible inhomogeneous mixing 
between the pristine gas and the AGB ejecta we adopt the 
following scheme. We assume that both the pristine gas and 
the ejecta are clumpy, and the number of these clumps or 
clouds is distributed in mass following a power law: 



dM 
dM 



(Al) 



(A2) 



where dMpr (dMoj) is the mass of pristine gas (ejecta) in- 
cluded in clouds with masses between M and M -\- dM. As 
a consequence, the total mass Mpr.tot of pristine gas, and 
Moj,tot of the ejecta, are given by 

A 



Mp 



Me- 



B 



{Ml-Z - Ml;-) 



(A3) 



(A4) 



where Mpr,u (Mcj,u) and Mpr,i (Mcj,i)are the lower and up- 
per extremes of the pristine gas (ejecta) distribution, and 
A and B are normalization constants. As in our program 
we form stars as a sequel of "bursts" occurring after every 
time interval At, Afpr.tot and A/cj,tot represent the amount 
of pristine gas and ejecta present during At, and are given 
by the simulation. As a consequence, taking A'/pr,u, Afoj,u, 
Mpr,i, Mej,i, a and /9 as free parameters, A and B are de- 
termined. Choosing an uniform deviate m^r in the range 
< JTipr < Afpr^tot we generate a random mass Afpr of pris- 
tine cloud from the distribution in eq. lAll and we proceed 
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simil arly to obtain a r andom value of the ejecta cloud Afej 
(e.g. iPress et al.lll992l ): 



A 



'2-/3 



pr,l 



l/(2-«) 



1/(2-/3) 



(A5) 



(A6) 



Assuming that these two clouds merge together, the 
resulting abundance Z in this mixed gas is given by 



(A7) 



where Zpr and Z^i represent the abundance of the element 
in the pristine gas and in the ejecta, respectively. 

After every At we chose a number of couples (Mpr, Afej) 
obtaining the .E-distribution of the stars forming at that 
time. 

As pointed out above, we have six free parameters that 
must be fixed. For illustration, we assume 



Afpr.i = 1 Mq . For the sake of simplicity, we further assume 
13 = a, Moj,u/Mej,i = Mpr,u/Mpr,i, and iVpr.tot = A^cj.tot (i.e. 
the number of pristine clouds is equal to the clouds of the 



ejecta), thus obtaining Moj,i/Mpr 



M 



ej jtot 



/Mpr 



which is given by the simulation at every At. 

In order to show the behaviour of the distribution 
for different values of a. Fig. lAll displays several his- 
tograms of the number of stars A'^([Na/Fe]) as a func- 
tion of [Na/Fe] forming in a mixture of two gases, one 
with [Na/Fe] =0, representing the pristine gas, and one with 
[Na/Fe] =0.82, representing the AGB ejecta. Two ratios be- 
tween the masses of these two gases are taken into ac- 
count: X = 1, characteristic of models like those in Paper 
I, and X = .05, typical of models as those suggested by 
IConrov fc Spe rgcl (2010), in which the original FG mass is 
the same observed today. 

We note that very steep (a = 3) and very flat (a = 0.1) 
mass functions give rise to similar distributions. For low val- 
ues of a the number of small clouds and massive clouds are 
similar, thus the probability of encounters between clouds 
with similar mass is elevated, and the final stellar distribu- 
tion peaks around the mean value which would be obtained 
in case of homogeneous mixing (vertical solid lines in Fig. 
lAlf) . In case of high values of a, the majority of the mass is 
distributed in small clouds of similar mass; merging between 
these clouds are the most probable, and once more the final 
distribution peaks at the mean value. For intermediate val- 
ues of Q, the larger clouds compensate their lower number 
with their larger mass. When they make the most probable 
merging with a small cloud, a large number of stars form, 
having a value of [Na/Fe] very skewed toward the value of 
the large cloud; for this reason A''([Na/Fe]) shows a clear 
bimodal structure with two peaks close to the edges of the 
[Na/Fe] range of values. 

The vertical dashed lines in Fig. lAll indicate the va lue 
[Na/Fe]=0.4 which is adopted by ICarretta et all (|2010l ) to 
separate the FG from the SG (see section (3)1 . The figure 
clearly shows that, for x = 0.05, inhomogeneous mixing can 
not populate the SG region ([Na/Fe] > 0.4) for any choice 
of a. This conclusion does not change even if we relax some 
of the other assumptions described above (and which play a 
minor role in shaping A''( [Na/Fe]). 
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Figure Al. DifTcrcnt distributions of A''([Na/Fe]) for different 
values ofa and x. The vertical solid lines indicate the mean value 
in case of homogeneous mixing (A''([Na/Fe]) would be shaped 
as a delta dirac function). The dashed lines separate the FG 
pop ulation ([Na/Fe]< .4) from the SG population ([Na/Fe] > 
0.4') l|Carretta et al.llSoiol 'l. For X = 0.05 the inhomogeneous mix- 
ing is notable to populate the SG region for any value of a. 



We investigate now the influence of a possible inho- 
mogeneous mixing models with e ~ 1, that is, in models 
in which the initial FG mass was much higher than to- 
day. Figure IA2I illustrates the fraction of SG stars with 
[Na/Fe]> 0.4 present in A''([Na/Fe]) as a function of x, for 
different values of a. It is clearly shown that, for large values 
of X, the fraction of new stars really appearing as SG stars 
to an observer oscillates in a range 0.6-0. 9, depending on a. 

From equation [T] (see text), the apparent ratio between 
the two populations is 



< 



Msg 
Mfg 



(A8) 



depending on the value of ^. Note that the maximum ratio 
is obtained for = 1, that is for the homogeneous mixing, 
when all the SG stars have [Na/Fe] > 0.4 (see the vertical 
solid lines in the left panels of Fig. lAll) . 
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Figure A2. Fraction f of SG stars with [Na/Fe]> 0.4 obtained 
by inhomogeneous mixing as function of x, and for different values 

of OL. 

Assuming ^ = 0.05 (see section [3]), it turns out from 
equation[T]that the apparent ratio between the masses of the 
two populations is a fraction ~ ^ of the true value (holding in 
the homogeneous case). As a consequence, models with inho- 
mogeneous mixing need the initial FG mass to be increased 
by a factor 1.1-1.7 with respect to homogeneous models to 
obtain the same amount of recognizable Msg. 
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